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ABSTRACT 

Laboratory searches for the detection of gravitational waves have focused on 
the detection of burst signals emitted during a supernova explosion, but have not 
resulted in any confirmed detections. An alternative approach has been to search for 
continuous wave (CW) gravitational radiation from the Crab pulsar. In this paper, 
we examine the possibility of detecting CW gravitational radiation from pulsars and 
show that nearby millisecond pulsars are generally much better candidates. We show 
that the minimum strain h c ~ 10~ 26 that can be detected by tuning an antenna to the 
frequency of the millisecond pulsar PSR 1957+20, with presently available detector 
technology, is orders of magnitude better than what has been accomplished so far by 
observing the Crab pulsar, and within an order of magnitude of the maximum strain 
that may be produced by it. In addition, we point out that there is likely to be a 
population of rapidly rotating neutron stars (not necessarily radio pulsars) in the solar 
neighborhood whose spindown evolution is driven by gravitational radiation. We argue 
that the projected sensitivity of modern resonant detectors is sufficient to detect the 
subset of this population that lies within 0.1 kpc of the sun. 



Subject headings: gravitation - radiation mechanisms: nonthermal - stars: neutron 
pulsars: general 



1. INTRODUCTION 



There are two types of signals of gravitational radiation that are expected to be most 
readily detectable from astrophysical sources: burst signals of short duration, from sources 
such as a stellar merger or a nonspherical core collapse associated with a supernova event; and 
continuous wave (CW) signals from sources such as short period, compact binary star systems or 
rapidly rotating, nonaxisymmetric (or precessing) compact stars. Most experimental searches for 
gravitational radiation have focused on the detection of burst signals. To date, there have been no 
confirmed detections of this type of gravitational radiation, although the terminal phase of the 
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coalescence of neutron-star binaries appears to be a promising source for future ground-based laser 
interferometers such as LIGO (Abramovici et al. 1992; Cutler et al. 1992; Finn & Chernoff 1993). 
The detection of CW radiation from binary star systems by a space-based interferometer has 
been proposed (Faller & Bender 1984; Evans, Iben, & Smarr 1987; Faller et al. 1989). Currently, 
though, the only experimental effort to detect CW gravitational radiation that is underway has 
been pioneered by the Tokyo group which has searched for CW emission from the Crab Pulsar 
(Tsubono 1991). We argue here that nearby millisecond pulsars are likely to be stronger sources of 
CW radiation than the Crab pulsar and therefore warrant the attention of experimental searches 
for gravitational radiation. 

An axisymmetric object which is rotating about its minor axis will not emit gravitational 
radiation because it has no time varying quadrupole moment. Therefore a pulsar (or any rotating 
neutron star, for that matter) must be nonaxisymmetric and/or precessing in order for it to 
radiate (Ferrari & Ruffini 1969; Zimmerman 1978, 1980; Shapiro & Teukolsky 1983; Barone et 
al. 1988). Several mechanisms for the production of nonaxisymmetric deformations in pulsars 
have been suggested (Ipser 1971), including asymmetric crystallization of the crusts (Ruderman 
1969; Ferrari & Ruffini 1969), pressure and magnetic stress anisotropies (Ostriker & Gunn 1969, 
Ruderman 1970), and rotationally induced instabilities (Imamura, Friedman, & Durisen 1985). 
Misalignments between the symmetry and spin axes of pulsars might occur as the result of 
electromagnetic torques due to magnetic dipole radiation, corequakes (Pines and Shaham 1974), 
or encounters with neighboring stars. 

Precessing, nearby millisecond pulsars have recently been put forth by de Araujo et al. 
(1994) as good candidates for the detection of gravitational radiation by upcoming interferometric 
detectors. They suggest, however, that the signals seen by these detectors due to radiation from 
wobbling pulsars may be burst signals, not CW signals. This is because the damping time-scale 
for the wobble angle due to the emission of gravitational radiation is expected to be on the order 
of seconds whereas the observation time needed to observe CW sources is likely to be ~ 10 7 s. 
This situation does seem likely to present itself with the first generation of detectors because the 
source signals are expected to be near the limit of detectability. Thus, at least for some time, 
nonaxisymmetric deformations may be the only channel through which pulsars produce detectable 
CW radiation. 

The rate at which gravitational energy is radiated from a nonaxisymmetric object that 
is rotating about its minor axis with angular velocity uj is (Ferrari & Ruffini 1969; Shapiro & 
Teukolsky 1983), 

£ G K = -ff/ 3 W. (1) 

This expression has been derived in the quadrupole approximation for nearly-Newtonian sources 
assuming that the object has principal moments of inertia I±, I2 and I3, respectively, about its 
three principal axes a ^ b > c fixed in the body frame (Landau & Lifshitz 1962; Misner, Thorne, 
& Wheeler 1973); and e = (a — b) /(ab) 1 / 2 is the ellipticity in the equatorial plane. Here, and 
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throughout this paper, dots denote differentiation with respect to time. Several authors have 
proposed that the Crab pulsar is the best candidate for detecting CW radiation (e.g., Zimmerman 
& Szedenits 1978; Tsubono 1991) because it has the largest spin down energy flux density 
Erot/{^r 2 ) of all pulsars; here r is the distance to the pulsar and E rot = I^ujuj (~ 10 38 erg s _1 
for the Crab pulsar) is the rotational energy loss rate. We note, however, that for a neutron star 
of given ellipticity e, the loss rate due to gravitational radiation strongly favors stars that are 
spinning more rapidly (equation [1]). Also, if the ellipticity is caused by the rapid rotation of the 
star, one might expect the value of e to be higher in stars that are rotating more rapidly. With 
this in mind, we have suggested (Barker et al. 1994; see also Schutz 1995) that nearby millisecond 
pulsars (Backer et al. 1982; Fruchter, Stinebring, & Taylor 1988; Johnston et al. 1993) may be 
stronger sources of CW radiation than the Crab pulsar (see §2). 

One favorable aspect of CW radiation is that a resonant detector can be tuned to the 
frequency of the emission. In order to tune a cylindrical bar detector so that its primary 
quadrupole mode of oscillation resonates with the frequency of radiation from the pulsar ujq = 2uj 
(see §2), the bar must have a length L « ttc s /ljq, where c s « 5.2 x 10 5 cm s _1 is the speed of sound 
in a prototypical aluminum alloy bar (see §3). Thus the length of the bar must be L ~ 1.3P ms 
m, where P ms = W 3 (2it/lu) is the period of the pulsar in milliseconds. This is prohibitively 
long (L ~ 43 m) for the Crab pulsar which has a period P ms = 33. Consequently, the Tokyo 
group (Tsubuno 1991) has used two short crossed bars instead of a single bar. Because better 
sensitivities are achievable using a single bar, millisecond pulsars appear to be better candidates 
from a detector design standpoint as well (see §3). 



2. THE CASE FOR MILLISECOND PULSARS 



In the quadrupole approximation, the two polarizations of the gravitational strain amplitude 
h that will be received by an observer located at a distance r along the rotation axis of a source 
are (Misner, Thorne, & Wheeler 1973; Kochanek et al. 1990; Rasio k Shapiro 1992) 

G 2- 

Here i xx ,i yy , and i xy are cartesian components of the reduced quadrupole moment, 

iij = hj - ^SijTrl, (3) 

where / is the inertia tensor and Tr/ = (I\ + I2 + h). If the source is a rotating, slightly 
nonaxisymmetric (a ~ 6, I\ 7^/2) star, 

Ixx - lyy = 4:L0 2 I 3 eCOs(2iOt) 
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I xy = 2iv 2 I 3 esm(2Lot) (4) 



(Shapiro & Teukolsky 1983) and 



G uj 2 

h+ = 4—r — l3ecos(2a;t) 
c 4 r 

G to 2 

h x = ^—hesin(2ut). (5) 

Therefore, the "characteristic" gravitational strain amplitude h c produced by a star that is 
spinning down due to the emission of gravitational radiation is 

G lo 2 

h c = 4^ — he (6a) 
c 4 r 

4.6 x 10 38 eros- i; wr' 1 ; 

where, in the second expression, we have inserted the value of I^e prescribed by equation (1) and 
the frequency at which the radiation will be detected is uo Q = 2oj. 



2.1. Expected Strains from Pulsars 



When a value of I3 = 10 45 g cm 2 is adopted (a reasonable estimate for a 1.4Mq neutron star 
of radius 10 km), equation (6a) can be written in the form 

^ = A.2 x 10^(PLn pc r\ (7) 

where r^ pc is the distance to the star in kiloparsecs. The right-hand side of equation (7) is given 
entirely in terms of observable pulsar parameters and assumes its largest value (h c /e ~ 10 -18 ) for 
nearby millisecond pulsars like PSR 1957+20 and PSR 0437-47 (see Figure 1). Unfortunately, one 
cannot use this expression to predict a particular strain from a given pulsar without making some 
assumption about the degree to which the pulsar possesses a nonaxisymmetric structure. 

Neutron star ellipticities could conceivably be as large as the breaking strain of their crust, 
which is thought to be ^ 5 x 10~ 4 (Shapiro & Teukolsky 1983). For any specific pulsar, however, 
an upper limit to e can be directly inferred from a measurement of P if one assumes that the 
spindown of the pulsar is attributed solely to the emission of gravitational radiation (Press & 
Thorne 1972). In this situation, E rot = I^uouo can be equated to Eqr as given in equation (1) and 
the inferred upper limit to the ellipticity is 

e GR = 6.0 (PLPf 2 , (8) 
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where, again, we have set I3 = 10 45 g cm 2 . Note that, as in equation (7), the right-hand side of 
this expression is given entirely in terms of observable pulsar parameters. Equation (8) serves 
only as an upper bound to the ellipticity because other mechanisms are likely to be responsible 
for some portion of the observed P in pulsars. Indeed, in most normal pulsars P is thought to be 
determined primarily by energy losses due to magnetic dipole radiation (Ostriker Sz Gunn 1969; 
Pacini 1967) rather than by losses due to gravitational radiation, so the true ellipticities of most 
pulsars will be less (perhaps orders of magnitude less) than the value egr given by expression (8). 
Note that the limit on e set by this expression for millisecond pulsars is very small; for example, 
for PSR 0437-47 and PSR 1957+20, e GK = 3x 10~ 8 and 1.6 x 10~ 9 , respectively (see Figure 1 and 
Table 1). 

It proves instructive to construct a plot of the two "observables" h c / e vs. egr for a large 
number of pulsars. Figure 1 includes (as of January, 1995) all of the (~ 470) pulsars in the 
Caltech Pulsar Database for which P ms ,P, and rk pc (the distance obtained from the dispersion 
measurement) were available. It should be noted, first, that because they tend to have very 
small values of P as well as small values of P ms , millisecond pulsars tend to cluster in the upper 
left-hand corner of this figure. By contrast, normal pulsars tend to cluster in the lower right-hand 
region of the diagram. The lines drawn in Figure 1, with slopes of —1, identify the locus of points 
for sources with identical values of h c rk pc as measured from Earth under the (unlikely) assumption 
that the spindown of all pulsars is attributed solely to the emission of gravitational radiation, i.e., 
under the assumption e = egr- For example, by this assumption, pulsars above and to the right 
of the solid line would have h c > lO -26 ^^ and the Crab and Vela pulsars, in particular, would 
exhibit the largest values of h c (~ 10~ 24 ). 

When confronted with a diagram like Figure 1, one is tempted to assume that the Crab 
and Vela pulsars will be the brightest CW sources in the sky. However, all that one can safely 
conclude is that, because these objects exhibit large spindown rates (which are understood to 
be due to magnetic dipole radiation, not gravitational radiation) the "observable" eqr places 
only a very loose constraint on the underlying structural ellipticity of these two neutron stars. 
One is permitted by eqr to adopt ellipticity values as large as few x 10~ 4 — and, hence, strains 
approaching 10~ 24 — without conflicting with measured values of P. However, in light of what 
we know about millisecond pulsars, it seems unreasonable to assume that these neutron stars have 
structural ellipticities of this magnitude. In particular, under the premise that nonaxisymmetric 
structure in neutron stars is likely to be rotationally enhanced (see Imamura, Friedman, &: Durisen 
1985), one would not expect e to be larger in (the relatively slowly rotating) normal pulsars than 
it is in millisecond pulsars. Correspondingly, it would be difficult to imagine that the Crab and 
Vela pulsars have ellipticities larger than e ~ 10~ 8 — 10~ 9 (the limit set by egr in millisecond 
pulsars) and, hence, that they produce strains larger than h c ~ 10~ 29 — 10~ 30 . 

In millisecond pulsars, it is conceivable that nonaxisymmetric distortions are induced by 
the observed, relatively large rotational energies of these objects and that measurable levels of 
CW radiation result. At first glance, PSR 0437-47 would appear to be the best candidate for 
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the detection of CW radiation because it is the millisecond pulsar with the largest value of the 
product (h c /e)eGR = 2.6 x 10~ 26 . However, taking the conservative view that e is nowhere larger 
than the value of egr set by PSR 1957+20 (i.e., e < 1.6 x 10~ 9 ), the pulsar exhibiting the largest 
h c is PSR 1957+20 itself (h c = 1.7 x 10~ 27 ). We consider this to be a much more reasonable 
estimate of likely strains coming from the best candidates for the detection of CW radiation than 
the estimates set earlier based on the observed properties of the Crab and Vela pulsars. 



In models that assume that the spindown of pulsars is due entirely to magnetic dipole 
radiation, it has been appreciated for quite some time that pulsar evolutionary paths can be drawn 
in an "observables" diagram similar to our Figure 1 if one adopts a particular function prescribing 
the strength of pulsar magnetic fields as a function of time (cf., Chanmugam 1992). Similarly, an 
evolutionary path in the h c /e vs. €gr diagram can be plotted for any individual pulsar if one 
assumes that the pulsar's measured P is entirely due to gravitational radiation and one adopts 
a particular function prescribing the pulsar's ellipticity as a function of time. To make such an 
evolutionary discussion relevant to the broad class of pulsars, rather than to one pulsar at a time, 
we should first remove any distance dependence from the diagram. By analogy with traditional 
discussions of stellar evolution, we define an "absolute" characteristic strain H c for any pulsar to 
be the characteristic strain it would exhibit if it were located a distance r^ vc = 1 from the sun. 
Figure 2 re-displays all the pulsar data from Figure 1 on this "absolute" characteristic strain scale. 

First, note that if P ms = 0.5 is associated with the maximum rotation rate (breakup velocity) 
of any neutron star, then the dotted horizontal line drawn at log 1 g (H c /e) = —16.78 in Figure 2 
represents the largest achievable absolute characteristic strain of any neutron star, even at birth. 
(The true demarcation line may, in fact, be somewhat higher or somewhat lower than this; its 
correct location cannot be established until our understanding of the equation of state of neutron 
star matter improves.) The simplest function to choose for eoRit) is cgr equals a constant. Then 
all evolutionary trajectories are vertical in Figure 2 and, because P is positive, all trajectories are 
directed downward. The rate at which a neutron star evolves along its vertical path in Figure 2 is 
prescribed by equation (8). Specifically, we can write 



2.2. Variation of Characteristic Strain with Age 




(9) 



where A = 1.9 x 10 5 s~ 3 / 2 . Then, assuming for simplicity that the pulsar period initially is much 
less than its period at time t, the function P(t) becomes 




(10) 
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Combining this with equation (7), the rate of evolution along a (vertical) trajectory in Figure 2 is 
prescribed by the following expression: 

^ = (4.0 x KT 19 * 1 / 2 ) e G l R t- 1 ' 2 . (11) 

According to equation (11), a population of neutron stars that have the same age but varying 
cllipticity will trace out a straight line that has a slope of -1 in Figure 2. More specifically, the 
solid and dashed line in Figure 2 locate isochrones of t = 10 s and 10 10 years, respectively. It 
is interesting to note that none of the individual pulsars plotted in Figure 2 lie below the 10 10 
year, or "Hubble line," isochrone. Somewhat surprisingly, isochrones in Figure 2 exhibit the same 
slope (-1) as lines of constant absolute characteristic strain (see Figure 1). Hence, a population 
of neutron stars of the same age will exhibit the same absolute characteristic strain across the 
entire population. For reference, the solid and dashed lines in Figure 2 correspond to strains of 
H c = 7.lx 10~ 27 and 7.1 x 1CT 28 , respectively. 



2.3. A Population of Loud CW Sources 



According to models that invoke magnetic dipole radiation to explain the spindown of pulsars, 
a pulsar's magnetic field strength is related to P and P by the expression B = 3.2 x 10 19 G 
s -i/2 (pp^j (Chanmugam 1992). By this relation, the extremely small value of P that has been 
measured for the millisecond pulsar PSR 1957 + 20 can be understood only if this pulsar possesses 
a magnetic field that is no stronger than 1.7 x 10 8 G. It seems unlikely that this field strength is 
unique among the population of rapidly rotating neutron stars (not all of which are radio pulsars) 
that resides in the solar neighborhood. Among this sub-population, gravitational radiation will 
dominate over magnetic dipole radiation as the principal energy loss mechanism throughout each 
star's life (<; 10 10 years) as long as the stars have equatorial ellipticities egr ^ 6 x 10~ 9 . (More 
generally, the criterion works out to be ecu <; 0.20-B 2 2 , where B\2 is the field strength in units of 
10 12 Gauss.) 

Given that the tightest constraint that currently has been placed on the ellipticity of neutron 
stars is e <J 5 x 10~ 4 , set by the breaking strain of the crust, it does not seem unreasonable to 
suggest that most neutron stars with B <J 10 8 Gauss also have structural ellipticities ecu <; 10 -8 . 
(As a point of reference, according to Cook (1973) the the gravitational field of the Earth exhibits 
an effective ellipticity ~ 10~ 6 .) We submit, therefore, that in the solar neighborhood there exists 
a family of rapidly rotating neutron stars whose spindown is driven by gravitational radiation. 
According to the preceeding discussion, even after 10 10 years of evolution we are guaranteed that 
this family of objects will exhibit an absolute strain ^ 7 x 10~ 28 . It would seem, then, that CW 
detectors designed to reach this strain level are virtually certain to detect this population of stars. 



-8- 



3. EXPERIMENTAL FEASIBILITY 

Experiments having the sensitivity necessary to directly detect gravitational radiation are 
difficult to design and conduct. We consider here a necessary condition for the feasibility of any 
new experiments designed to improve significantly the sensitivity for detection of CW sources. 
The specification of sufficient conditions is a much broader issue requiring discussion beyond the 
scope of this paper. 

It is well-understood that there are two fundamental, and predictable, noise sources that 
will necessarily limit the sensitivity of any gravitational wave detector: thermo- mechanical force 
noise and sensor noise. The first of these arises from random forces applied to the antenna due 
to its coupling to a finite-temperature heat bath; the second results from random fluctuations at 
the output of the electro-mechanical devices that are used to sense any gravitationally induced, 
time-varying distortions of the antenna. For detection of burst sources, both types of noise are 
important (Solomonson et al. 1993). 

The situation is different if one conducts a targeted search for CW emission from a known 
pulsar. Then information from radio observations can be used to great advantage. For example, 
from the measured pulsar rotation frequency to, one knows to what frequency uj = 2uj the 
antenna should be tuned; and because P is extremely small, one expects to receive a steady flux 
of gravitational radiation at the prescribed frequency for a very long time (P/P ~ 10 9 — 10 10 yr 
for millisecond pulsars). The continuous interaction of a resonantly tuned antenna with a steady, 
periodic gravitational force allows the mechanical amplitude of the antenna to steadily grow with 
time. Coherent averaging of the motion, accounting for the observer's and source's doppler shifts, 
can be used to narrow the bandwidth of the measurement to such an extent that the sensor noise, 
under reasonable conditions, can become completely negligible. Thus it should be possible to 
avoid sensor noise altogether within a narrow range of frequencies and thereby greatly reduce the 
minimum detectable strain amplitude h c compared to the amplitude that is detectable from burst 
sources. Such techniques have been employed in the search for CW radiation from the Crab pulsar 
(Tsubono 1991). 

In practice, then, the mechanical motion of the antenna will have two parts: the motion 
induced by the gravitational wave signal and the motion induced by the thermo-mechanical force 
noise. This remaining noise source is simply the random Langevin force F, or the mechanical 
equivalent of Johnson- Nyquist noise in electrical circuits. The relative amplitude of these two 
contributions to the signal can be calculated by starting from the equations of ordinary elastic 
theory and including the CW gravitational and Langevin forces as driving terms (Merkowitz & 
Johnson 1994). The elastic equation can be solved via an eigenfunction expansion which produces 
one harmonic oscillator equation for each eigenmode of the system. In this harmonic oscillator 
formulation, the CW gravitational force takes the form 

fGW = \^Leh + {t) = --flL e hcUo cos u Q t, (12) 
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where fi is the effective mass of the detector and L e is the effective length of the mode for a 
particular component of the gravitational strain tensor. Representing the effects of F by a random 
gravitational field h R , the Langevin force f F (t) takes the form 

f F (t) = ^L e h R (t). (13) 

Therefore the total force acting on the antenna can be written as: 

hot = ^fiL e (h R (t) - h c u 2 cosuj t). (14) 

The spectral density S R of any random force f R is given by 

Sr =< f* R f R >, (15) 

where f R is the Fourier transform of f R and an asterisk denotes the complex conjugate. Setting 
Ir = If, and realizing that h R = —co 2 h R , we deduce 

S F = (^L e u 2 ) 2 S h , (16) 

where 

S h ^< h* R h R > . (17) 

But the single-sided spectral density of the Langevin force near a resonant mode at frequency uo 
is fixed by simple thermodynamics and the fluctuation-dissipation theorem to be 

S F = 4kT(^), (18) 

where k is the Boltzmann constant, T is the physical temperature of the heat bath, and Q is called 
the mode quality factor (Q~ l is the dissipation coefficient). Hence, combining equations (16) and 
(18) we derive 

This expression gives the spectral density of the noise produced in the antenna by the Langevin 
force, normalized to the amplitude of the force on the antenna due to the gravitational wave itself. 

The square root of Sh is referred to as the "strain noise" of the antenna and can be used 
to define the background noise level below which a signal of astrophysical origin will remain 
undetectable. More specifically, if one integrates the signal over an observing time r and demands 
a 4cj confidence level for detection, the smallest reliably detectable h c will be 

(20) 



When developing a detector for CW gravitational radiation, expression (19) is important 
because it identifies design parameters that contribute most significantly to the amplitude of the 
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random noise near the resonant mode of the antenna. It is important to realize, in addition, that 
the factors of length, mass, and frequency are not independent. They are coupled by the geometry 
and composition of the antenna. For example, for the primary quadrupole mode of a "long" right 
circular cylinder, or bar, with length L and diameter d « L, the mode parameters become 

p = -pird 2 L, (21a) 
8 

4 Ar 

L e ~ -rL ~ — , (216) 

where p is the mass density of the antenna material and c s is the speed of sound in the material. 
Substituting these expressions into equation (19), the expression for strain noise becomes 

JS~ h = (J^LS) 1 ' 2 = 1.2 x 10~ 23 ( — ) ( — ! I — ! 11/ (22) 
V \cPp4Qj \ d J \0.05KJ { Q ' K J 




where we have used the density (p ~ 2.7 g cm~ 3 ) and sound velocity of aluminum alloy 
(c s ~ 5.2 x 10 5 cm s _1 ), which we regard as the most likely material for such an antenna. Note 
that neither uj nor L has a direct effect on the sensitivity of the bar antenna, but the source 
frequency does affect the length L of the bar (L = 2.1 m for PSR 1957+20, for example). 

With current technology, the chosen scale parameters are within the realm of feasibility. 
The specified diameter, for example, is at the limit of current commercial casting capability, and 
T ~ 0.05-fT has recently been achieved for the NAUTILUS bar detector under development by 
the Rome group (Astone 1991). We have assumed that a moderate improvement in the value 
of Q over presently available bar antennas can be accomplished. Such values of Q have been 
achieved by annealing small samples of commercially obtained aluminum alloy 5056 (Marsden 
1984). The most significant technical difficulty in achieving a strain noise as low as 10 -23 Hz -1 / 2 , 
as indicated by our parameterization of expression (22), is likely to be vibration isolation in such 
an ultra-low-temperature cryostat. 

An integration time ~ 4 months seems feasible for tracking a pulsar that is well timed 
from radio observations, and makes it reasonable to repeat the measurement a number of times 
for confirmation. Combining expressions (20) and (22), and setting r = 10 7 s, the smallest 
detectable h c becomes 1.5 x 10 -26 . By comparison, the Tokyo group has obtained an upper limit 
of h c ~ 2 x 10~ 21 for the Crab pulsar (Owa et al. 1988). Current experiments they are conducting 
are likely to lead to improvements of a factor 100 or so (K. Tsubono, private communication). 

A recently proposed, alternative to the cylindrical bar detector is a "spherical" detector 
having an antenna configured as a truncated icosahedron, that is, an antenna with the same 
geometric shape as the Cqq Buckyball molecule (Johnson & Merkowitz 1993; Merkowitz & Johnson 
1994). An antenna of this type with a diameter D exhibits 5 resonant quadrupole modes, each 
having mode parameters of the form: 

w «3.24^, (236) 
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H = p^D 3 , (23a) 
L e » 0.301D. (23c) 

Again, the size of the antenna is fixed by the sound velocity and the source frequency; for example, 
D = 2.2 m for PSR 1957+20. Substituting these expressions into equation (19) gives the following 
expression for the strain noise of a spherical antenna: 

The first equality appears to be little different from the bar case but, actually, because the 
diameter D for a sphere is larger than the diameter d for a bar, when both have the same resonant 
frequency the value of \J~Sh for the sphere is smaller. Additional advantages of a spherical detector 
are that it can be equally sensitive to a wave from any direction and it is capable of measuring 
the direction and polarization of the wave (Forward 1971; Wagoner & Paik 1976; Merkowitz & 
Johnson 1994). 

In the second equality of expression (24) we have used the coupling of the antenna's diameter 
to the frequency to eliminate D in favor of the pulsar period P ms . We have again assumed 
aluminum alloy as the material. As in the case of the bar detector, a significant technical difficulty 
will be maintaining high Q when casting or joining aluminum to make a sphere of diameter 2.2 m 
designed to resonate with twice the rotation frequency of PSR 1957+20. If the integration time is 
again chosen to be r = 10 7 s, the minimum detectable strain for a spherical detector tuned to PSR 
1957+20 is predicted to be h c = 8.0 x 10~ 27 . 



4. DISCUSSION 



Nearby millisecond pulsars are good candidates for the detection of CW gravitational 
radiation. Because of their close proximity and rapid rotation, they are capable of emitting 
radiation with larger gravitational strain amplitudes than pulsars that are farther away and/or 
have longer periods. 

The minimum strain h c ~ 10 -26 (1.5 x 10~ 26 for a bar detector, 8.0 x 10~ 27 for a spherical 
detector) that can be detected by designing an antenna tuned to the rotation frequency of the 
millisecond pulsar PSR 1957+20 and employing presently available resonant detector technology is 
orders of magnitude better than what has been accomplished so far by observing the Crab pulsar, 
and within an order of magnitude of the maximum strain that can be produced by PSR 1957+20 
as a result of rotationally induced nonaxisymmetric deformations. The design and operation of a 
resonant antenna that is tuned to the rotation frequency of PSR 1957+20 would, at the very least, 
place physically meaningful constraints on the nonaxisymmetric ellipticity of millisecond pulsars. 
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We have argued (§2.3) that there almost certainly is a population of rapidly rotating neutron 
stars (not necessarily radio pulsars) within the solar neighborhood whose spindown evolution is 
driven by gravitational radiation. Throughout their entire lifetime, these stars will radiate at 
an "absolute" strain H c ^ 7 x 1CT 28 . It is significant that the projected sensitivity of modern 
resonant detectors is sufficient to detect the subset of this population of stars that resides within 
0.1 kpc of the sun. 

Support through NASA grant NAGW 2447 and NSF grant PHY-9311731 is gratefully 
acknowledged. 
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Table 1: Data for Representative Pulsars 3, 



Pulsar 


p 


P 


^kpc 


h c /e 




Crab 


33.4 


4.21E-13 


2.49 


1.52E-21 


7.55E-4 


Vela 


89.3 


1.25E-13 


0.61 


8.70E-22 


1.80E-3 


PSR 0437-47 


5.76 


1.2 E-19 


0.14 


9.12E-19 


2.89E-8 


PSR 1957+20 


1.61 


1.68E-20 


1.53 


1.07E-18 


1.59E-9 



"From the Caltech Pulsar Database 
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Fig. 1.— The two "observables" h c /e = 4.2 x Kr 18 (P£ s r fcpc ) _1 versus e GR = 6.0 (P^ S P) 1/2 , for 
~ 470 pulsars from the Caltech Pulsar Database. Binary pulsars are denoted by asterisks. The 
solid and dashed lines identify the locus of points for sources with identical values of h c rk pc = 10~ 26 
and 10 -27 , respectively, under the assumption that gravitational radiation is solely responsible for 
the spindown of the pulsars, and therefore, that e = €gr- 
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Fig. 2. — The same data shown in Figure 1, except that all pulsars are now assumed to be located 
at a distance r^ vc = 1 from the sun. On this scale, the "absolute" characteristic gravitational strain 
amplitude is H c = 4.2 x 10 _18 P~^e. If the breakup velocity limits the period of a neutron star to 
Pms = 0.5, the dotted horizontal line marks the largest possible H c /e for any neutron star. The 
solid and dashed lines represent the isochrones of pulsars with ages of 10 8 and 10 10 (the Hubble 
line) yrs, respectively, under the assumptions that €gr is a constant function of time and e = eqr- 



